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ABSTRACT 

The dynamical friction and diffusion coefficients are derived for a massive binary that moves against 
a uniform background of stars. The random impulses exerted on the binary's center of mass by the 
field stars are greater than those exerted on a point particle due to inelastic scattering. The frictional 
force acting on the binary is less than that acting on a point particle due to randomization of the 
trajectories of field stars that pass near the binary. Both effects tend to increase the random motion of a 
binary compared with that of a point mass. If the maximum effective impact parameter for gravitational 
encounters is comparable to the radius of gravitational influence of the binary, its Brownian velocity can 
be increased by a modest factor (< 2) compared with that of a single particle. This condition is probably 
fulfilled in the case of binary supermassive black holes in galactic nuclei. 

1. INTRODUCTION 

Brownian motion, the irregular motion exhibited by a heavy particle immersed in a fluid of lighter particles, is usually 
modelled by the Langevin equation, 

d i = -Av + F(t). (1) 

Here v denotes the velocity of the heavy particle, — Av is the frictional force from the fluid, and F(t) is the fluctuating 
force that arises from collisions between the heavy particle and the fluid particles. The precise form of F(t) is typically not 
known but its statistical properties are well defined; for instance, the mean value of F is zero and its fluctuations generally 
occur on a time scale that is short compared to the characteristic time over which v varies. When these conditions arc 
met, the evolu tion of the single-particle distribution function / is described by the Fokker-Planck equation (e.g 



Kampen 1992); the dynamical friction coefficient that appears in that equation is just A, and the diffusion coefficient is 
related to the amp litude of F via the requirement that the steady-state motion be in energy equipartition with the fluid 



at temperature T ( Einstein 19551 ), 

This paper is concerned with tnc Brownian motion of a massive binary, consisting of two point particles whose center of 
mass moves irregularly due to random gravitational encounters with field stars. The Brownian motion of a binary differs 
from that of a single particle because the kinetic energy of a field star is not conserved following a close interaction with 
the binary. If the binar y's mass is m uch greater than that of a field star, most binary- field star interactions will extract 
energy from the binary (Heggie 1975), increasing the field star's kinetic energy and hence the recoil velocity of the binary's 
center of mass. This "super-elastic scattering" will give the binary a larger random velocity than expected for a point 
particle in energy equipartition with background stars, or 



M 



(2) 



where M and m/ are the mass of the heavy particle and of a field star respectively, and (v"j) is the mean square velocity 
of the field stars. 

A motivation for the work presented here is the likely existence of binary black holes. There are at least two environments 
where binary black holes are expected to form. In globular clusters, stars with initial masses exceeding ~ 20M© should 
comprise ~ 0.1% by number of the first generation of stars; these massive stars would leave b ehind black holes with masses 
M, ~ IOMq, which would rapidly fall to the cluster center and form black-hole binaries ( Portegies Zwart fe McMillan 
1999). Subsequent close enc ounters will eject most of the bla ck holes, with perhaps a single binary remaining at the 
cluster center after ~ 10 9 yr ( Portegies Zwart fc McMillan 2000). Bin ary black holes are also expected to form in galactic 
nuclei following galaxy mergers ( |Bcgclman, Blandford fc Rccs 1980| ) ; here the characteristic masses are M, ~ 10 8 M Q . 
The evolution of a binary supermassive black hole due to encounters with stars has been modelled by Makino (1997), 
Quinlan & Hernquist (1997) and others using A-body codes. Quinlan & Hernquist (1997) found that the wandering of 
the binary in their simulations was a factor 5 — 10 greater in amplitude than expected on the basis of equation (|^), and 
attributed this difference to inelastic scattering of stars off the binary. 

More speculatively, the massive halos of the Milky Way and other galaxies may contain large numbers of black holes. 
A population of stellar-mass, primor dial black holes ha s been proposed to explain the gravitational microlensing results 
toward the Large Magellanic Cloud ( Alcock et al. 1997). A small fraction of these black holes would be expected to form 
binaries through three-body interactions QNakamura, Sasaki fc Tanaka 1997). Alternatively, the halos of ga laxies like the 
Milky Way may consist of black holes of mass ~ 10 D Aifn ( Carr, Bond fc Arnctt 1984 ; Lacey & Ostrikcr 1985) which would 
sink toward the galactic center and form binary systems there ( Hut fc Rees 1992| ) 
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The contribution of inelastic scattering to the random motion of a massive binary can be estimated as follows. Define 
a to be the semi- major axis of the binary and M 12 = Mi + M 2 its mass, with M 12 3> rrif. The impact parameter p 
corresponding to a closest-approach distance r p of a field star to the binary is 



P 



1 



2GM 15 
V 2 r n 



IGMxi 



(3) 



with V the relative velocity at infinity. The rate of encounters with stars that interact strongly with the binary r p < a, 

is then ~ nf(irp 2 )a f w 2imfGM\ 2 aj a j with rif the number density of field stars. These stars are ejected with typical 

velocities v e j ~ y GM\ 2 /a : the binary orbital velocity (Hills & Fullerton 1980), yielding a rate of energy extraction from 
the binary of 



dE b 



dt 



with pf = rnfUf. We can also write 



1 

2 m/ " 



2nn f GM 12 a\ nG 2 M{ 2 p f 



dE b 



dt 



d (G{M 12 /2) 2 



dt 



2n 



(4) 



(5) 



(assuming M\ = M 2 ), whic h, together with equation 
constant rate (Heggie 1975): 

d_ 

dt 



implies the well-known result that a hard binary hardens at a 
8irGpf 



The hardening rate is usually expressed in terms of a dimensionless constant H as 



d_ 

dt 



H 



Gpj_ 



Three-body scattering experiments (Hills 1983; Mikkola & Valtoncn 1992; Quinlan 1996) give H 
(GM 12 /a ^> a 2 ), equal-mass, circular-orbit binary with M\ 2 3> m/. 
In terms of H, the energy extraction rate is given precisely by 



dE h 



dt 



d 
~dt 



GM 2 2 



G 2 M 2 2 p f 
8a f 



H. 



(6) 
(7) 

15 for a hard 
(8) 



Momentum conservation implies that almost all of the binding energy released during an encounter will go into kinetic 
energy of the field star. Thus 



SE h w SE„ 



1 



mf{8v^) 2 



1 /Mia. 
-m f I Sv 



V m f 



(9) 



with Sv the change in the center-of-mass velocity of the binary; the last relation follows from conservation of momentum. 
The rate of diffusion in velocity of the binary due to super-elastic scattering is then 



((Av) 



)S.E. 



2m/ 
m l2 



dE b 



dt 



H G 2 m sPf 
4 a s 



This diffusion rate may be compared with Chandrasekhar's expression for a point mass initially at rest, 

(N2)c=8(2l) i/^l nA , 



(10) 



(11) 



with In A = In (p m ax / Pmin) , the Coulomb logarithm ( Spitzer 1987 , eq. 2-12). Both diffusion coefficients scale in the same 
way with the parameters (mf,pf,af) that define the field star distribution, except for the weak variation implicit in the 
Coulomb logarithm. The ratio of the two coefficients, 



H 



32V27T In (pmax/Pmin) 



(12) 



therefore depends only on In A and H, reflecting the fact that the Chandrasekhar coefficient is dominated by distant 
encounters while the effects of super-elastic scattering are limited to close encounters. A more exact definition of In A, and 
appropriate choices for p max and p m in, will be presented below; here we note that p m in is of order GM\ 2 /a 2 = tq, the 
radius of gravitational influence of the massive binary. If the distribution of field stars around the binary falls off steeply 
over a radius comparable to tq , p max will be comparable to p m i n and In A will be small. Equation (03) therefore predicts 
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that super-elastic scattering will increase the amplitude of the binary's random velocity by a modest amount compared 
to the random velocity of a point particle with the same mass. 

There is a second way in which energy exchange between field stars and the binary will act to increase the magnitude 
of the binary's random velocity. Brownian motion represents a balance between dynamical friction and scattering. But 
field stars that interact strongly with the binary are ejected in nearly random directions, and this reduces the dynamical 
friction force that they exert on the binary. The velocity change experienced by a field star in a low-impact-parameter 



collision with a point-mass perturber is ~ — 2V, corresponding to a 180° change in its direction (cf. equation 28). When 
the point mass is replaced by a hard binary, the field star is ejected in a nearly random direction and its mean velocity 
change (averaged over many encounters with different phases and orientations of the binary) is therefore ~ — V in a 
direction parallel to V. The drag force exerted on the massive object is proportional to the mean velocity change of the 
field stars and hence the contribution to the frictional force from close encounters is only ~ 1/2 as great in the case of a 
binary as in the case of a point mass. 

At first sight, this reduction in the frictional force might appear to be significant. In the standard treatment (e.g. 



Spitzer 1987, eq. 2-46), dynamical friction arises entirely from field stars with velocities less than that of the massive 
object. If a massive binary were moving slowly, all of these field stars would have small relative velocities as well and 
would interact strongly with it. However it turns out (§2.1) that the standard treatment is incorrect in this regard and 
the reduction in dynamical friction is correspondingly more modest (§4.1). 

The approach adopted in the present paper consists of calculating the dynamical friction and diffusion coefficients for a 
massive binary under exactly the same approximations adopted by Chandrasekhar (1942) in his analysis of gravitational 
encounters with a point mass. In the case of a binary, this calculation requires the numerical integration of a large number 
of binary- field star interactions, i.e. scattering experiments (§4). The dynamical friction and diffusion coefficients can then 
be computed by taking appropriate averages over the velocity changes experienced by the field stars in these experiments. 
Comparing the coefficients computed for the binary to those computed by Chandrasekhar for a point mass gives the 
increase in the expected magnitude of the binary's random velocity compared to that of a point mass with M = M\%- 
The restriction in the present paper to a massive binary, M\2 ~S> mj, permits a number of simplications; for instance, the 
mass of the field star can be neglected in the scattering experiments, and the values of the dynamical friction and diffusion 
coefficients need only be computed in the limit v <^ <jf appropriate to a massive body that is near energy equipartition 
with much lighter stars. 

It became clear after this investigation was begun that the Chandrasekhar coefficients as they are usually presented 
could not be directly compared with the results of the scattering experiments. Starting with Chandraskehar (1943a), the 
dependence of the Coulomb logarithm on the velocity of the field star has generally been i gnored when integrating over 



the field star velocity distribution (e.g. Elosenbluth, MacDonald & Judd 1957; Henon 1973). This approximation is valid 
in the case of a test star whose velocity is comparable to that of the field stars but not when v <C 07, which is the case 
of interest here. Thus in §2 the Chandrasekhar coefficients are re-derived without the usual approximations. One of the 
interesting results of this re-derivation is that the dynamical friction force on a slowly-moving body is found to come 
primarily from field stars with Vf > v. In the standard treatment, all of the frictional force comes from field stars moving 
more slowly than the test star. 

An important preliminary step is to relate the equilibrium velocity distribution of a test mass (single star or binary) to 
its dynamical friction and diffusion coefficients, via the Fokker-Planck equation. This is done in §2. Applications of the 
results to various physical systems are discussed in §5. 

2. FOKKER-PLANCK EQUATION 

We wish to compute the equilibrium velocity distribution f(v) of a massive binary that interacts via gravitational 
encounters with individual stars; v is the velocity of the binary's center of mass. We assume that the time scale over 
which this equilibrium is reached is much shorter than the time scale over which the orbital elements of the binary 
change; this assumption is discussed in more detail in §5.2. We also ignore a possible dependence of / on the direction 
of v, appropriate if the orientation of the binary relative to v is rapidly changing and if the field stars are distributed 
isotropically in velocity. The dependence of the dynamical friction force on the orientation of the binary is discussed briefly 
in §4.1. Because of the assumed large mass ratio, 3> rrif, even close encounters will produce only sm all deflections in 
the trajectory of the binary. Hence we may describe the evolution of / via the Fokker-Planck equation ( phandrasekhar 



1943b). In a steady state, this equation implies 



1 d 



= -N(v)(Av) + -—[N(v)((Av) 2 )} (13) 

where 



4irv 2 f(v)dv = N(v)dv. (14) 

The coefficients (Aw) and ((Aw) 2 ) are defined in the usual way as sums, over a unit interval of time, of Aw and (Aw) 2 due 
to encounters with field stars. 

Define Auii and Awj_ to be changes in w in directions parallel and perpendicular to v. We have 

r 2 i 1 / 2 

Av= (v + Avn) + (Aw ± ) 2 -v. (15) 
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Expanding to second order in the small quantities Avh and Av± and taking means, the coefficients in equation (|l3| ) may 
be expressed in terms of the usual dynamical friction and diffusion coefficients as 



\2\ 



(Av) = (Av l{ ) + ~ {(AV v ±J ((At,) 2 ) = ((A W|| ) 2 ). (16) 



Substituting equations ( |l4| ) and (|T^) into equation (|T^) yields^] 



= / 



{Av ] \) + ^[{Avl)-2{Avl)) 



Equation ( jl7j ) may also be derived from equation (2-71) of Spitzer (1987) after replacing dE in that equation by vdv and 
expressing [AE) and (AE 2 ) in terms of (Av\\), (Av 2 ) and (Avj_). 

In a steady state, we expect that v will be of order \Jrrif /M12 <C 1 times the velocity dispersion a f of the field stars. 
Hence it is appropriate to expand the dynamical friction and diffusion coefficients about v = 0: 

(Av l{ ) = -Av + Bv 3 . . . , (18a) 
(Avf) = C + Dv 2 . . . , (18b) 
(Avl) =2(E + Fv 2 )... (18c) 

Substituting these expressions into equation (|l7]) gives 

= (-Av 2 + Bv 4 + E + Fv 2 -C- 2Dv 2 ) f - I (C + Dv 2 ) ^. (19) 
v ' 2 v ' av 

Setting v — gives C = E, i.e. the diffusion rate of a stationary particle must be independent of direction. Cancelling 
additional terms, we find to lowest order in v 

= (A + 2D-F)f+~?£K (20) 
2 av v 

The standard Chandrasekhar coefficients for a point particle of mass M interacting with stars of mass to/ (e.g. Spitzer 
1987, eqs. 2-52 - 2-54) give, after expanding about v — 0: 

4V2^G 2 m 2 f n f liiA ( M\ 3 A c . , 

3 aj \ m f J 10 cry 



8V2^ G 2 m 2 n f \nA 3 C c , 

Cc = ^ ^ ' ^"lO^' 

1 C c 



Thus 



E C =C C , F c = (21c) 



A c + 2D C -F C =A C -^ (22a) 



7 

iV2n G 2 Mm f n f In A 
~3 ^ 



(22b) 



« A c (22c) 

where the last relation is valid when M 3> rrif. In the case of the more general coefficients describing the motion of a 
massive binary, A will still be of order Mi 2 /to/ times D and F and we may likewise approximate (A + 2D — F) by A. 
Thus 

OvAf + ^l (23) 
2 av v 

which has solution 



f(v) = he~ v2 ' 2 °\ (24a) 
2A 

In equation (rLTj) and below, the parentheses in expressions like ((Any) 2 ) are dropped 



(24b, 
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Substituting Aq for A and Cq for C gives the well-known result for the equipartition velocity dispersion of a point mass, 

It will be useful in what follows to express the coefficients A and C for a massive binary in terms of their values for a 
point mass with M = M12, i.e. in terms of the ratios 

fl lS A (26 ) 

(The Chandrasekhar coefficients Ac and Cc will be defined slightly differently below.) The velocity dispersion of the 
binary is then 



73 \ 1/2 / \ 1/2 

#2^ ( m f ^ 



RiJ \M 



'12 



Thus the problem of determining the equilibrium velocity distribution of the binary's center-of-mass motion is reduced 
to determining the ratios of (A«||) and (Av 2 ), the dynamical friction and diffusion coefficients for the binary, to the 

Chandrasekhar coefficients for a point mass in the low-v limit. We note that R2 can equally well be defined as (C + 
2E)/(C C + 2E C ), i.e. in terms of the diffusion coefficient (Av 2 ) = (Av?,) + (Air[), in the limit v -> 0. 



3. CHANDRASEKHAR COEFFICIENTS FOR A POINT MASS 



The Chandrasekhar coefficients for a point mass in the form given above (eqs. 21a-21c) embody two approximations: 
the argument of In A was assumed constant when integrating over the distribution of field-star velocities; and the non- 
dominant terms, which are of order 1/lnA or 1/A times the dominant terms, were ignored. Both approximations are 
reasonable when the velocity of the test star is not too different from the typical velocity of the field stars, and when 
the integration over impact parameters extends to distances much greater than GM/a 2 ,. However it is impossible to 
reproduce these approximations precisely when carrying out the scattering experiments discussed below; hence in order 
to compute the ratios (eqs. |2^) between the coefficients for the binary and for a point mass, we must first rederive the 
Chandrasekhar coefficients without the usual approximations. Of particular interest here is the velocity dependence of 
the Coulomb logarithm. Chandrasekhar (1943a) and White (1949) noted that this dependence implies that some of the 
dynamical friction force comes from field stars moving more rapidly than the test mass; in the standard treatment, where 
In A is removed from the integral over field star velocities, all of the friction is produced by field stars with Vf < v. An 
exact analysis in the case of interest here, v <C 07, has apparently never been carried out. We will find in fact that almost 
all of the frictional force comes from field stars with Vf > v in the low- velocity limit. We also derive a precise expression 
for the Coulomb logarithm in this limit. 

3.1. Dynamical Friction 

The velocity change of a test particle of mass M in one encounter with a field particle of mass mj « M is flSpitzci 



1987, eq. 2-19) 



Av\\=-2V r ^-- ^ (28) 

11 M l+p 2 /pl 

where V is the relative velocity at infinity, p is the impact parameter, and pg = GM/V 2 . Multiplying by 2irpnfVdp, with 
n/ the number density of field stars, and integrating over p gives 

N = - M ^ h( 1+y L/FS, (29) 

the rate of change of V\\ due to encounters with field stars of velocity V. 
The dynamical friction coefficient is (Spitzer 1987, eq. 2-28) 



(A«,|) = j f f (v f ) (At;,,) ^"^dv f (30a) 

= -2nG 2 Mm f n f J f f ( yj )?LZ?LZ. ]n (l + §ff^) dv,, (30b) 

with v the velocity of the (massive) test particle, // the distribution of field-star velocities (normalized to unit total 
number) and V = v — ; v has been assumed parallel to the rr-axis. Following Chandrasekhar (1943a), we can represent 
the velocity-space volume element in terms of vj and V; then writing 
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gives 

poo rv + Vf 

(Av l{ ) = -27r 2 G 2 Mm f n f v- 2 / dvj v f f f (v f ) / dV [ 1 + — — >- J In ( 1 + ) (32) 



\V-Vf\ 




where ff is henceforth assumed isotropic 



In the standard approximation (e.g. Rosenbluth, MacDonald & Judd 1957), the velocity dependence of the logarithmic 
term is ignored, and one writes 

l..(l + f^)-2,„A, (33) 
The dynamical friction coefficient then becomes 

(A«||) = -4ir(4irG 2 Mm f n f ) [ dv f (^-Y f f (v f )Hx (v,v f , p max ) (34) 

Jo \ v / 

with 



if v < Vf. 

Neither p max nor p m i n are well-defined; p m i n is typically set to some multiple of GM/a 2 (e.g. Spitzer 1987, eq. 2-14). 
Equation (|3^) reproduces the well-known result that only field stars with Vf < v contribute to the frictional force. 

When the velocity of the test particle is sufficiently low, as in the case of a massive particle near equipartition with 
lighter field particles, the logarithm in equation (|32|) will be close to zero for all stars with Vf < v. Hence it is unclear 
whether removing the logarithm from the integrand is a reasonable approximation, or what the proper definition of In A 
(particularly p m in) in the final expression should be. The frictional force in the low-w limit must either be much less than 
implied by equation (|34|), or else most of the frictional force must come from field stars with Vf > v. 

In fact the latter is the case, as we now show. Returning the logarithmic term to the integrand in equation (EmJ) gives 



H 1 (v,v f , Pmax ) = — dV [1 + ——L In 1+^— . (36) 




This function is plotted in Figure la. For low test-star velocities, v < \J GM/p max , Hi deviates strongly from a step 

function; it reaches a peak at Vf » y/ GM/p max and its peak amplitude varies as ~ v 3 due to the velocity dependence of 
the logarithm. Figure lb plots the (normalized) integrand of the dynamical friction integral, equation (|34|), in the limits 
of low and high v/at. The field-star velocity distribution was assumed to be a Maxwellian: 

//(»/) 'j^n^- <W 

For low vjcxf, the integrand of equation (p3) may be shown to vary approximately as 



x f e~ x V 2 
1 + i? 2 4 



(38) 



where 



_v f p max _ GM 

x f = —> R = Pf = — §-• 39 

of pf aj 

This function peaks at Xf w i?" 1 / 2 . In the case of high v/af, on the other hand, the integrand varies as 

x)e- x V 2 (40) 

which peaks at Xf = 1. Thus for R = p ma x/pf ~ 1, the dynamical friction force comes from field stars with Vf ~ (jf, 
whether the test star is moving much faster or much slower than the typical field star! For R > 1, the frictional force in 

the low-v limit comes from field stars with Vf < Of/R 1 / 2 ; thus for any v < af/R 1 ^ 2 , i.e. for M > Rrrif, most of the stars 
producing the frictional force will be moving faster than the test star. The appropriate choice for R is discussed in §5 
where it is argued that R > 1 in many cases of interest. 

The complete expression for the dynamical friction coefficient is 

(A«||) = -16^/irGmfPfnfF(v/af,p max /pf), (41a) 

F(x,R)=x~ 2 dyy 2 e- y2 H 1 (x,V2y,R), (41b) 
Jo 

r x+z / t 2 — r 2 \ 

(41c) 



1 f x+z ( x 2 - z 2 \ 

Hi(x,z,R) = — dw[l + ^—)\n(l + R 2 w 4 ) 

8zJ\x-z\ V w ) 
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where x = v/cxf and R and pj are defined above. The function F(v,p max ) is plotted in Figure 2. The dependence 
of (Auii) on v is seen to be a function of R. By contrast, when In A is assumed fixed, the w-dependence of (Auy) is 
independent of p max /p m i n (Spitzer 1987, eq. 2-52; Fig. 2). 



The low-f limit of the frictional force may be derived by careful manipulat ion of equations (41a - 41c). A more 
transparent route, useful in what follows, is to first change variables in equation (30a). Q Representing the veL 
volume clement in terms of V and A, where A is the direction cosine between V and v 



ocity-space 



A = 



(ve x - Vf) ■ ve x _ v - v fiX 
Vv V [ 



gives 



(Avn 



poo pi \ / 2 jr4 

-2irG 2 Mm f n f x 2tt jf dVV 2 dX— f f (v f ) In (l + 



-4Tr 2 G 2 Mm 



( p 2 V 4 



dWff(vf). 



Again substituting equation ( |37| ) for ff and integrating over A, 

2v/2^G 2 Mto 



(A«,|)=- 



7 



V 2 V 2 



vV 



vV 



vV 



— 7T cosh — — sinh — =- 



7 



7 



7 



Taking the limit v <C 07, and writing by analogy with equation (18a) 

(Av\\) = -A c v + B c v 3 

we find 

2^2^ G 2 Mm f n f 



Ar 



Be 



7 



2?r G 2 Mm f n f 



15 



dz e- z ln(l + 4RV 



poo 

/ dz e~ 2 (5- 2z)ln(l + 4i? ; 



Thus Hooke's law ((Awy) cx — v) is recovered for u <C cry. 
Of primary interest here is the leading coefficient, which may be written 

Ac = ^G*Mrn f n fG{R)> 



i r°° 

G(R) = - / dz e- z ln(l + 4i? : 
A reasonable approximation to G(R), valid for R > 1, is 



G(R) « In ^l + 2^«lnWl+^^ 



(42) 

(43a) 
(43b) 



(44) 

(45) 

(46a) 
(46b) 



(47a) 
(47b) 



(48) 



(Figure 3). Comparing equations fl47a| ) a nd (M8 |) with the standard dynamical friction coefficient for a point mass in the 
low-test-particle-velocity limit, equation fl21a|), we find 



]uA = ku/l 



Prnax 



or 



In VT + 2i? 2 = lnWl 



GM 
V2cr 2 ' 



^Prnax° f 



max u f 

G 2 M 2 



(49) 
(50) 



J I am grateful to M. Milosavljevic for pointing this out. 
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Spitzer (1987, eq. 2-14) advocates 



G(m f + M) 



(51) 



in the case that the test and field stars have similar masses and velocities. When M » rrif, the denominator in Spitzer's 
expression should be reduced by a factor <~ 2 to account for the lower relative velocities, giving a value for p m in that is 
not too different from the value derived here. This agreement is not fortuitous: it is a result of the fact that the field stars 
responsible for producing the frictional force have roughly the same velocity distribution whether the test star is moving 
rapidly or slowly. 

3.2. Diffusion 

The squared velocity changes of the test particle in one encounter with a field star are (cf. Spitzer 1987, equations 2-18 
and 2-19): 



Multiplying by 2wpnfVdp and integrating over p as before, 



2_,T,2 m / 1 V_ A1 r2™f P 2 /P 2 



(Av ± Y = AV 2 



M 2 (l+p 2 /p 2 ) 2 ' 



(Avn) 2 = 



(Avl) 



2 _ 



AirG 2 nfm 2 
V 

4irG 2 nfm 2 
V 



Pmaxl Pp 
1 +Pmax/PoJ ' 



In II 



Pi 



Pmax/Pa 
l +P 2 maxlPl 



(52) 

(53a) 
(53b) 



the rates of change due to encounters with field stars whose velocity at infinity is V. 

The expressions just given refer to velocity changes with respect to the direction of the initial relative velocity vector 
V. We need to transform to a frame in which the test particle has velocity v. Following Spitzer (1987, eqs. 2-24, 2-25), 
but retaining the contributions from the non-dominant terms, gives in the fixed frame 



1 



(AviAvj) = (e, • ei)(ej • ei)(A«||) 2 + - [(e< • e^e, • e' 2 ) + (e, • e^)(e, • e' 3 )} (Av ± ) 2 



V 2 



or 



(A Vll r - -(A V± r 



V, 



+ MA V± r 



v 



(Av 1 Av 1 ) = -^(Av«) 2 + 7 : (At, 



(Av 2 Av 2 ) + (Av 3 Av 3 ) 



V 2 
V 2 + V 2 
V 2 



V 2 , 

/ iV 2 + V 2 ^ 



(Avj_) 2 



where the e! axis is oriented parallel to v. The final integrations over field-star velocities are then 

fV+V f 



(Aujj) = — / dv f Vff f (v f ) / dV V (AuiAwi), 

V Jo J\v-v f \ 

(Av 2 ± ) = — dv f Vfff(v f ) / dVV [(Av 2 Av 2 ) + (Av 3 Av 3 }\ . 

v JO J\v-v f \ 



Writing 



we find as before 



V x = 



V 2 + v 2 - v 2 



f 



2v 



Vy' + V 2 = 1- V 2 , 



(A«f) = ^-{A-nG 2 m)n f )v dv f (^f f f (v f )H 2 (v,v f ,p, 

8ir f°° /Vf\ 2 

(A«i) = —{A-KG 2 m)n f )v dv f (^J f f (v f )H 3 (v,v f ,p, 

3 r v + v f 
H 2 (v,v f ,p max ) = - — / dV 

OVf Ji v _ v I 



Av 2 



2 2 

v — vj 



v 2 V 4 
In I 1 + Pmax 
1 + G 2 M 2 



(54a) 
(54b) 

(55a) 
(55b) 

(56a) 
(56b) 

(57) 

(58a) 
(58b) 

(58c) 



Merritt 




H 3 {v, Vf,p max ) = / dV^ 




P 2 ta ^ 4 /G 2 M 2 
l+P 2 max VyG*Mi> 

2" 



y2 

pL^ 4 /g 2 m 2 

l+^a^ 4 /G 2 M 2 



V 2 ^ 4 

In 1 + ^ ma3: 
1 G 2 M 2 



(58d) 



In the standard approximation, In A is assumed to be a constant and the non-dominant terms are ignored, yielding 

l) r ~- r r . 

h^ = { , , r:\ . f (59) 



lnA(^f) z if <;>«,, 
lnA(^) i£v<v f . 



In A (3 
21nA 



if v > «/, 

if V < Uy. 



(60) 



We are again primarily interested in the low-u limits. Returning to equations (56a) and (56b), expressing the integral 
in the same (V, A) variables as in §3.1, and expanding to second order in v/af gives 



(Avf) = C C + D c v 2 . . . , 
(Avl) = 2{E C + F c v 2 ) 



(61a) 
(61b) 



C c 


8 , G 2 m 2 f n f f°° 

= W2k f dz e^lnfl- 


- AR 2 z 2 ) , 




(62a) 


D c 


2 . — G 2 m 2 f n f f°° 

V2^ / / dz e~ z (5 
15 g) Jo V 


+ 2z-2z 2 ) ln(l 


+ AR 2 z 2 ) , 


(62b) 


Ec 


= C C , 






(62c) 


F c 


2 , G 2 m 2 f n f r°° 

V2^ I 1 / dz e~ z (5 
15 a) J V 


-6z + z 2 ) In (l -+ 


AR 2 z 2 ) . 


(62d) 



The leading terms are 



Cr — Ec — 



/2vr G 2 m 2 f n 



G(R), 



(63) 



with G(R) again given by equation (47b). Comparing equation ( J63| ) with equation (47a), we see that the ratios between 
the three diffusion coefficients at low v are exactly the same as in the standard treatment with In A assumed constant. 
Thus we predict the same equilibrium velocity dispersion for a point mass as in equation (p5|), a 2 — Cc/2Ac — (mf /M)a 2 . 

In what follows, the expressions (47a) and ([33]) will be used to define Ac and Cc when computing the ratios R\ and R 2 
of the dynamical friction and diffusion coefficients of the binary to those of a point mass (eqs. ||^) . 

4. SCATTERING EXPERIMENTS 

The dynamical friction and diffusion coefficients for a massive binary were computed by carrying out a large number 
of scattering experiments, in the manner of Hills (1983) and Quinlan (1996). In the limit of large mass ratio Mu/rrif, 
the fi eld star may be treated as a massless particle moving in the potential of the two black holes (Mikkola & Valtoncn 
1992), which follow unperturbed Keplerian trajectories about their fixed center of mass with semimajor axis a. From the 
changes in the field star's velocity during the encounter we can infer the corresponding changes in the binary's center-of- 
mass motion. Let the initial velocity of the field particle be Vf = v fo^x and its final velocity (v' f e x + v' f e y + v' f e z ). 

The velocity changes (Sv x , dv y , 5v z ) of the binary's center of mass following a single collision are 



Sv x = 


m f I i \ m f x 
M J Vf * V ^ = M 12 SVf *> 


(64a) 


5Vy = 


mf v f > = mf Sv f , 
Mia J y Mia ly 


(64b) 


Sv z = 


m f . ' m f X . 
M 12 VJz ~ Mia Vfz ' 


(64c) 
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The dynamical friction and diffusion coefficients for the binary are obtained by taking the appropriate averages, over an 
ensemble of trajectories, of the Sv^s and their squares, as discussed in more detail below. 

Orbits were integrated using the routine D0P853 of Hairer, N0rsett & Wanner (1991), an 8(6)th order embedded Runge- 
Kutta integrator. The routine automatically adjusts the integration time step to keep the fractional error per step below 
some level TOL, which was set to 10~ 9 . The forces from the black holes were not softened. Each field star was assumed 
to begin at a position (x,y, z) = (oo,p, 0) and was advanced from r = oo to r — 50a along a Keplerian orbit about a 
point mass M%2- The integrations were terminated when the star had moved a distance from the binary that was at 
least 100 times its initial distance with positive energy, or when the number of integration steps exceeded 10 6 . The latter 
condition was met only for the integrations with the lowest initial velocities, and even then, for much fewer than 1% of the 
integrations. These stars were on orbits that were weakly bound to the binary and that needed to make many revolutions 
before being expelled. The incomplete integrations were not included when computing the diffusion coefficients below. 

The binary was given unit total mass, Mi = M 2 = 1/2 and its center of mass was fixed at the origin. The binary's orbit 
was assumed to be circular in all of the integrations described below. The orientation of the binary's orbital plane, and 
its initial phase, were chosen randomly for each integration. The angle £ between the normal to the binary orbital plane 
and the initial relative velocity was stored for later use, since velocity changes in the field star are expected to depend 
systematically on £. However most of the results presented below are based on averages over £ and over the phase of the 
binary, appropriate if the orientation of the binary relative to v is changing rapidly due to encounters. 

Initial velocities for the scattering experiments were assigned one of the discrete values K x Vt>m = K x yj GM\ 2 /a 
with K = (0.01,0.02,0.03,0.05,0.1,0.2,0.3,0.5,1,2,3,5,10). 10 6 orbits were integrated for each value of K, for a to- 
tal of 1.3 x 10 7 integrations. The impact parameters p for each V were chosen from one of 20 intervals correspond- 
ing to ranges in scaled pericenter distance r p /a of [0,0.001], [0.001,0.003], [0.003,0.01], [0.01,0.02], [0.02,0.03], [0.03,0.05], 
[0.05,0.1], [0.1,0.2], [0.2,0.3], [0.3, 0.5], [0.5, 1], [1,2], [2, 3], [3, 5], [5, 10], [10, 20], [20, 30], [30, 50], [50,100] and [100,200]. Im- 
pact parameters were chosen randomly and uniformly in p 2 within each interval. Unless otherwise indicated, distances 
and velocities given below are in program units of a and (GM12/ 'a) 1 / 2 respectively. 

A number of checks of the numerical integrations were carried out. The field-star velocity changes were used to 
recompute the binary hardening rates calculated by Mikkola & Valtonen (1992), Hills (1983,1992) and Quinlan (1996). 
The agreement was excellent. The relation: 

4v 



Hi(V) = Q2 ^ 2n |V(A«,|) + ((A«n) 2 + (A V J 2 )] , (65) 



was also checked; here Hi is the velocity-dependent hardening rate defined by Quinlan (1996), and the diffusion coefficients 
on the right hand side refer to the field stars. A weaker check (since it does not depend on the details of binary-field 
star interactions) consisted of verifying that the diffusion coefficients (Adm) and (Air^) for the binary tended to those of 
a point mass for large p ma x ■ 

4.1. Dynamical Friction 

For a binary initially at rest, the acceleration of its center of mass due to encounters with field stars of initial velocity 
V is 

(<MV,iw)) = ^^J^ d P P ^x(V.p) = -27m/ C^j vj^ ma * dp p foJ~ x (V, P ) (66) 



where Svf (V, p) is the mean change in the cc-component of the field star velocity, averaged over many encounters with 
the binary at fixed (V,p). We are interested in the ratio between the numerically-computed coefficients and the exact 
expression for a point mass (equation ^) . This ratio is 

Di(V,jw) = - 7-7-1 , „ 2 7I2Y7^T72 / d P P Sv fx ( 67a ) 



^(l+P 2 ma JP 2 o)G 2 M 2 J 

V \ 3 [ Pmax /a , (p\ 2 ffrv* 



(!) & 



2 In (1 +p 2 max /pl) \Vbbx) J \aJ \V h 

with V£ m = GM 12 /a and Po = GM 12 /V 2 . 

Figure 4 shows values of D\, averaged over all orientations of the binary, for several values of V. Generally D\ < 1 
due to randomization by the binary of the velocities of test stars that pass within a distance r p < a. However there is a 
complicated dependence of D-y on p max and V. For V ^ Vbim collisions at small impact parameter (p < a) result in the 
field star passing through the binary with almost no change in velocity on average, and hence D\ — > 0. As V is decreased, 
the velocity changes in the field star become greater on average and D\ is determined by the distribution of angles over 
which field stars are ejected from the binary. Figure 5 shows the distribution of field star velocity changes as a function 
of p for V — Q.hVbin- At large impact parameters, p 3> p C rit, where 



p crit 2GM 12 . . 

~ w (68) 
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the velocity changes are distributed in the way expected for interaction with a point mass (Fig. 5a; eq. g8|). When 
P ~ Pcriu the distribution of Svf's broadens due to the change in velocity (both magnitude and direction) resulting from 
interaction with the binary; the broadened distributions extend to positive Svf's, corresponding to stars which are ejected 
with high velocity in the same direction as their initial motion - producing an "anti-frictional" force on the binary (Fig. 
5c). At still smaller impact parameters, p < p cr it, stars are ejected with a typical angle ~ 7r/3 with respect to their initial 
velocity vector, compared to ~ tt for interactions with a point mass; thus | is reduced and with it the contribution 
to the frictional force (Fig. 5d). 

The dynamical friction coefficient for the binary is (cf. equation |30a| ) 

(A«||) = f f f (v f ) ( V -^-) (Sv x )dvt. (69) 



We are primarily interested in the low-v limit. Repeating the analysis leading to equation (47a) we find, for a Maxwellian 
.//• 



i /2 , f°° jv \ / y x ! 



This quantity was computed after scaling the velocity changes obtained from the scattering experiments, given assumed 
values for the two dimensionless parameters 

Ae^, S = — , (71) 

Pf °f 

with pf = GMi2/<J 2 t\ the parameter S measures the hardness of the binary. Figure 6 shows R\ = A(R, S)/Ac, the 



dynamical friction coefficient for the binary normalized by the point-mass coefficient (equation 47a). Values of p m ax ^ a 
are not physically interesting, hence R and S are restricted to R > l/S 2 . The greatest reduction in the frictional force 
occurs at small R and S. Small R corresponds to low p m ax, hence to a small contribution of distant encounters to (Ai>ii). 
Small S implies a wide binary (in comparison with p max , say) and hence a larger fraction of stars that interact strongly 
with the binary. In the case R « 1, i.e. p max ~ Pf, we see that the reduction in the frictional force is ~ 50% for the 
widest binaries, Vbin ~ 07, falling to ~ 20% for Vu n — 2a f and ~ 5% for Vun = 5°7- 

Results presented so far were averages over all orientations of the binary's orbital plane with respect to the initial velocity 
of the field star, and the direction of the binary's center-of-mass motion. In reality, the binary will have a particular 
orientation with respect to its center-of-mass motion, at least on time scales short compared to the time required for the 
orientation of the binary or its direction of motion to change. The typical velocity change of an interacting field star would 
be expected to depend systematically on the angle VP between the normal to the binary orbital plane and V. Figure 7 
illustrates the dependence Di(^,p max ) for the case V = 1.0 X y/GMw/a. The dynamical friction force is greatest when 
the binary orbital plane is perpendicular to the relative velocity vector, although the dependence on $ is mild. In the 
case of low impact parameters, p max < a, and VP ~ 0, the frictional force is close to zero since the field stars pass through 
the middle of the binary at high enough velocity that almost no deflection takes place. 

4.2. Diffusion 

We may similarly define dimensionless ratios D2 and D3 between the diffusion coefficients (Sv 2 ) and (5v 2 +5v 2 ) computed 



from the binary scattering experiments, and their equivalents for a point-mass scatterer, equations ( p3a|) and (53b) 




(72a) 



(72b) 

Figures 8 and 9 plot D2 and D3 as functions of p max and V. The coefficient D2 describes mean square changes in the 
binary's velocity parallel to V; these are non-dominant, i.e. bounded in their integration over impact parameter at large 
Pmax- Hence D2 does not tend toward unity at large p ma x- At small p m ax, the expected change in Vf for scattering off 



a point mass is ~ — 2V (equation 52), whereas a typical change after interaction with the binary is ~ —0.5y/GMi2/a (cf. 
Fig. 5). These are comparable when V « 0.25Vbj„; thus for V < 0.25Vfejn one expects D2 > 1 for all p ma x, while for 
V > 0.25Vbin one expects D% < 1 for all p ma x, in good agreement with Figure 8. For sufficiently large V, D 2 drops to 
zero at small p ma x since most field stars pass through the binary without significant change in velocity. 

-D3 describes the mean square cha nge i n the binary's velocity perpendicular to V, which diverges as lupmax in the case 



of a point-mass perturber (equation 53b ); thus D% — ► 1 in the limit of large p m ax- For small p, the change in Vf for a 
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point-mass scatterer goes to zero (cf. equation 52), while Svf 2 remains finite in the case of interactions with the binary; 
hence D 3 diverges at small p max . 

The diffusion coefficients for the binary are 



V\ ( V x 2 



2E =^r d &)(^) 2 '- v *' 2 '' lSv2 ' +Sv ->- <73b) 

Figure 10 plots the ratio R2 — (C + 2E)/(Cc + 2Ec) between (Aw 2 ) for the binary and for a point-mass perturber in 
the low-w limit (cf. §2). The dashed line in that figure is 

1 -1 H (74) 

32v^7rlnVl + 2i? 2 

the expected ratio in the limit of an infinitely hard binary according to the approximate treatment of §1, equation (p"2|); 
In A in that ex pression was replaced by ln-y/l + 2R 2 (eq. [19J) and the hardening rate H was set to 15, appropriate for 



a hard binary (Hills 1983; Quinlan 1996). The computed curves may be seen to tend toward the analytic expression as 
S = Vbin/cff increases, i.e. as the binary becomes harder. 

Finally, Figure 11 gives the ratio a/ac = (R2/ Ri) 1 ^ 2 , the factor by which the equilibrium velocity dispersion of the 
binary is increased relative to that of a point mass with M — M.\%. Based on equation (|74|), this ratio is expected to be 




°- 18 , (75) 
lnVl + 2i? 2 V ' 

in the limit of a hard binary; this estimate ignores the additional contribution to the random motion from the reduction 
in dynamical friction discussed above. Equation (f75j) is plotted in Figure 11. At fixed Pmax, Figure 11 implies that ojoc 
will increase with increasing hardness of the binary, until reaching a limiting value similar to equation ([75|). Thus the 
wandering of the binary should increase slightly in amplitude as its semi-major axis shrinks. 

The results of this section may be summarized as follows. The velocity-dependent dynamical friction coefficient D\(V) 
for a binary (Figure 4), with V the velocity of the field star at infinity, is typically less than that for a point particle due to 
randomization of the direction in which field stars are ejected during close encounters (Figure 5). The velocity-dependent 
diffusion coefficient D^(V) (Figure 9) is typically greater than that for a point particle due to the gain of energy which 
a field star experiences during a close encounter. After averaging over field-star velocities, the dynamical friction and 
diffusion coefficients for the binary, (Aum) and (Au 2 ), are respectively less (in absolute value) and greater than their values 
for a point particle (Figure 6 and 10). Both effects act to increase the random velocity of a binary in equilibrium with field 
stars compared to the equilibrium random velocity of a point particle of the same mass (Figure 11). However the increase 
is modest unless R = Pmax/Pf = Pmax&'j /GM12 is less than about one, in other words, unless the maximum effective 
impact parameter for gravitational encounters p max in Chandrasekhar's theory is of order the radius of gravitational 
influence of the binary GMn / a 2 . 

5. DISCUSSION 
5.1. Estimating p max /pf 

The results presented above, in particular equation J75| ) and Figure 11, demonstrate that the random motion of a 
massive binary will be dominated by distant (elastic) encounters and therefore essentially the same as the Brownian 
motion of a point particle, unless the ratio 

2 

p _ Pmax PmaxVf 

R = — = ^Lj (?6) 

is smaller than about one. The motivation for this result was presented already in the Introduction: the net effect of 
distant (elastic) encounters is proportional to In A w In Vl + 2i? 2 (cf . equation [l8|) , hence R must be small for close 
(inelastic) encounters to dominate. 

If the massive binary is located in a constant density region of radius ~ r c at the center of a stellar system with a 
steeply falling density profile, then p max ~ r c (Appendix A). 

To estimate pf, we note that pf ~ rc, the radius of gravitational influence of a body of mass Mi 2 embedded in a stellar 
system with velocity dispersion 07. In most applications of Chandrasekhar's theory, pj is smaller than p max by orders of 
magnitude and R S> 1. For instance, scaling to parameters appropriate to stars in a globular cluster, 

GM Grnt , „ / 07 \ 2 / mf\ 

p f = w — w 4.3 x 10 J pc i r — f- , (77) 

; 0) (t 2 ^ \10kms / \ M &) 
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much smaller than a core radius of ~ 1 pc. However in the situations of interest here, pf is larger, and R smaller, by 
roughly a factor M 12 /mf ^> 1 (ignoring changes in ay). In the case of a binary consisting of two 10M Q black holes in a 
globular cluster, we have 



Pf w 10 3 pc 



°7 



/ Mi 



,10 km a -1 / V 2OM 0, 

In the case of a supermassive black hole binary at the center of a galaxy, we can use the M, — a relation 

4.7 

M i2 = M, ~ 1.30 x 10 B M Q ' 



°7 



200 km s" 1 



(Merritt fc Ferrarese 2001) to write 



Pf Rj 10 pc 



°7 



200 km s" 



2.7 



(78) 



(79) 



(80) 



Consider first a binary black hole at the center of a globular cluster. A core radius of r c ~ 1 pc is still much greater 
than pf w 10~ 3 pc, implying _R ^> 1 and no appreciabl e enhancement of the Brownian motion. However in the "post- 
core-collapse" globular clusters ( Djorgovski fc King 1986| ) the stellar density follows a steep power law all the way into the 



observable center and r c is much smaller than 1 pc. While p max is not well defined in this case, arguments like those in the 
Appendix suggest that the effective p ma x would be small, hence R < 1 and a substantial enhancement in the Brownian 
motion is predicted (Figure 11). 

Next consid er the case of a ~ 10 8 M^ black-hole binary in a galactic nucleus. Li ke globular clusters, some galaxy nuclei 
have "cores" (Lauer et al. 1995), actually weak power laws in the space density (Merritt fc Fridman 1996). (The cores 



of globular clusters could also be weak power laws.) The "break radius" r u in these galaxies jFabcr et al. 1997 ) plays 
roughly the role of a core radius; Poon & Merritt (2000) note that r& « a few x tq in two galaxies where both n, and 



M. can be accurately measured, hence p m ax <; Pf and R > 1. In the "power-law" galaxies (Ferrarese et al. 1994; Lauci 
et al. 1995| ), which have no detectable core, R would be smaller, R < 1. The enhancement of the Brownian motion could 
be substantial for a massive binary at the center of a power-law galaxies, perhaps as much as a factor of two in a (Figure 

The concept of a maximum effective impact parameter p ma x on which these estimates depend is a loose one, both as 
used here and in Chandrasekhar's original theory. Refining these estimates of p ma x will require A-body simulations. 

5.2. Equilibration time scales 

Here we test the assumption made above (§2) that the center-of-mass motion of the binary reaches a statistical steady 
state in a time short compared to the time over which the orbital elements of the binary change. In fact this is not always 
the case. The time scale for change of the binary's semi-major axis (cf. equation is 



tha 



dt 



Gp f aH 



Approach to statistical equilibrium of the Brownian motion occurs in a time (cf. equation 18c) 



The ratio is 



with S 2 = V 2 



tha 



I ha 



V 




v 2 


~ 1 ~ 1 _ 3 CT / 


(A«,|) 




(Aw 2 ) 


~ A ~ A c 4V2^G 2 Afi 2 p/lnA 



4V2^ GM 12 



^relax 



3H 



In A 



4V2^S 2 In A 3 . 35 2 In y/l + R 2 



'I 



3H 



H 



(81) 



(82) 



(83) 



/(7 2 = GM 12 /(T 2 a. This ratio must be large if the binary's semi-major axis is to remain nearly fixed 



during the time required for the Brownian motion to be established. For large enough R, i.e. large enough p r . 
guaranteed. However the case of most interest here is R » 1 (§5.1), for which 



tharde 



relax 



h ' 



this is 



(84) 



For a very hard binary, H w 15 (|Quinlan 1996| ) and t& eca y/t re i a x ~ 0.07V^ 2 in /cr 2 , which exceeds unity for S = V\,inlo~f ^ 4. 
Thus the condition is satisfied for sufficiently hard binaries. For a soft binary, S < 1, we can use the expressions in 
Quinlan (1966, eqs. 16-18), to write H w 4AS 2 , giving tdecay/tielax ~ 0.23. Thus in the case of a soft binary, the binary 
separation will change appreciably during the time that the Brownian motion is being established. However in the soft 
regime, the Brownian velocity dispersion depends only very weakly on S (cf. Figure 11), so the conclusions presented 
above are not significantly affected. 
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5.3. Consequences of Brownian motion 
The Brownian velocity dispersion of the binary is predicted to be (equation ^7|) 



i?i Mi 



(85) 



We can convert this into a mean square displacement by assuming that the binary moves in a harmonic potential well 
produced by a uniform density core of field stars. The virial theorem 



= -(F-r) 



(86) 



as applied to the binary's center-of-mass motion relates the mean square velocity of the binary to an average of the force 
acting on it. In a spherical constant-density core, the force (neglecting encounters) is — (AirGpf /3)r and 



Writing r 2 = Qa 2 f /^Gp s ( |King 1966| ) yields 



(v 2 ) = -vG Ps (r 2 ). 



R 2 m-f 



(87) 



r 2 Ri M u 

This treatment (cf. |Bahcall fc Wolf 1976 ) is approximate in that it ignores the effect of the motion of the binary on the 
background stars and assumes that gravitational encounters are uncoupled from the quasi-periodic motion in the core. 
For a binary black hole at the center of a globular cluster, the rms Brownian velocity is 



and the wandering radius is 



Vrms ~ 3.87 km s 1 



0.22 pc 



R2 

Ri 



1/2 /on \ 1/2 
' / 20m - ■ 



V M 



12 



/ 



10 km 



R 2 \ 1/2 [2Qm f \ 1/2 I r, 



^Rx) V M i2 / V 1 P c , 
In the case of a binary supermassive black hole in a galactic nucleus, the Brownian velocity is only 

V 2 / ™ \ l/ 2 



0.030 km s 



and the wandering radius is 



0.0087 pc ( 



1/2 



RiJ \200kms" 1 



200 km s" 



-2.35 



-1.35 



100 pc 



(89) 



(90) 



(91) 



(92) 



the M. — a relation ( |79| ) has again been used to relate M 12 to cry. The Brownian motion of a supermassive black hole is 
of course very small, even with the enhancement due to inelastic scattering derived here. The predicted wandering radius 
is smaller than the likely separation of the black holes in a binary syst em; in fact it is of order the separation at which 
gravitational radiation would lead to a rapid coalescence (Merritt 2000) 



Measureme nt of the velocity of the Milky Way black hole has recently become feasible (Reid et al. 1999; Backer & 



Sramek 1999). The upper limit of 
3 x 10 6 M o 



20 km s 1 is a factor ~ 200 greater than the expected Brownian velocity given its 



mass of • 

The wandering a mplitude of a binary supermassive black hole in a galactic nucleus has r ecently been discussed by a 
number of authors (Makino 1997; Quinlan & Hernquist 1997; Merritt 2000; Gould & Rix 2000). The problem is interestin g 
because a stationary black-hole binary will quickly eject all of the stars with pericenters less than ~ 2a (e.g. Zier 2000); 
once this occurs, the density of stars around the binary drops and the decay stalls. Wandering is a possible way to increase 
the number of stars that a bina ry can intera ct with, thus allowing the binary to decay to the point that gravitational 
radiation coalescence can occur (Peters 1964). 

Quinlan & Hernqust (1997), in a series of computer simulations using a hybrid iV-body code, noticed a wandering of a 
massive binary with an amplitude more than five times greater than expected on the basis of an equation like (|92|). The 
core radius that they cite, ~ 0.04, is comparable to tq in their simulations, so R « r c /ra ~ 1- Figure 11 suggests that 
the increase in the amplitude of the wandering would probably be much less than the claimed factor of 5 — 10. 

Makino (1997) carried out TV-body merger simulations similar to those of Quinlan & Hernquist (1997) but using a more 
conservative, direct-summation code. Makino's Figure 7 shows a wandering amplitude that scales as ~ TV -1 / 2 ; the rms 
velocity of the binary in that plot is comparable, or perhaps ~ 50% larger, than expected for a point mass, consistent 
with the results obtained here. 
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Quinlan & Hernquist (1997) used a spectrum of masses for the field particles, spanning a range of ~ 10 3 in mt, and 
adopted the lowest mass when computing the expected amplitude of the binary's Brownian motion. Mass segregation 
might have brought some of the more massive particles into the nucleus over the course of their simulations. Quinlan 
& Hernquist also used a basis-function-expansion code for computing the forces between field stars. Such codes do not 
conserve momentum and this too may have contributed to the random motion of the binary. A repeat of the Quinlan & 
Hernquist simulations using a different iV-body code could help to clarify the discrepancy between their results and those 
of Makino. 

If supermassive black hole binaries wander as little as implied by equations ([75|) or (|9|), it is difficult to see how 
they could int eract with enough stars to allow them to achieve gravitational-radiation coalescence in a reasonable time 
( |Merritt 200C| ). Unless some other mechanism is effective at removing energy from these binaries, one might expect them 



to persist for the lifetime of a galaxy at separations a < 0.1 pc, or at least until a subsequent merger event brings another 



supermassive black hole into the nucleus (Valtoncn 199£) 



I thank Milos Milosavljevic for advice about evaluating some of the integrals and for numerous discussions about the 
ideas presented here. A conversation with Douglas Heggie helped to motivate this paper. This work was supported by 
NSF grants AST 96-17088 and 00-71099 and by NASA grants NAG5-6037 and NAG5-9046. 

APPENDIX 



Maoz (1993, eq. 4.4) gives an implicit expresson for p max appropriate for an object at the center of a spherically- 
symmetric matter distribution p(r). Maoz's formula, 



p(o)y Pm ,„ r 



dr. 



replaces the Coulomb logarithm in expressions like equation ( 21a ). If we set 

„2\ -7/2 



p{r) = p(0) 1 + - 



with r c the core radius of the field stars, Maoz's formula gives 



In A 



sinh 



if 7 = 1, 



Setting p m i n w pf Rj GM\ija\ (eq. p0|) and equating In A with In \/l + 2R 2 (eq. 49), we find 



p f r G 



(Al) 



(A2) 



(A3) 



(A4) 



for both values of 7. 
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Fig. A2. — The dimensionless dynamical friction coefficient F(v,p max ) (equation 4ft) for p m ax/pf = {0.3, f , 3, f 0, 30, f 00, 300, f 000}. 
Dashed curve is the classical expression with In A = 4; vertical dotted line indicates the position of the maximum of this function. 
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Fig. A3. — The function G(R) (equation 



47b). The dashed line is the approximation of equation (|48|), G(R) w In vT-j- 2R' 2 . 
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Fig. A4. — Dynamical friction reduction factor Di(V, Pmax) (eq. B7b). D\ is denned as the mean change in field-star velocities after 
interaction with the binary, expressed as a fraction of the value expected tor interactions with a point-mass scatterer of the same total mass as 
the binary; p m ax is the maximum impact parameter and a is the binary semi-major axis. Different curves correspond to relative velocities at 
infinity of V/V^i n = 0.1, 0.2, 0.5, 1, 2, 10; the lowest (highest) V produces the highest (lowest) D± at large pmax- Points are averages computed 
from the numerical integrations; curves are spline fits. Some of the curves have been extended to large p ma x/a> using the analytic expression 
for a point-mass scatterer. 
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Fig. A5. — Distribution of field star velocity changes for scattering experiments with V/Vt,i n = 0.5. Each plot corresponds to 5 X 10 4 
scattering experiments within some range of impact parameters [pi,p2] in units of a. (a) [6,10] (b) [2,4] (c) [0.6,1] (d) [0.4,0.6]. Solid lines in 
(a) and (b) are the distributions corresponding to scattering off a point-mass perturber. In (c) and (d), the mean of this distribution (which 
is very narrow) is indicated by the arrows. The magnitude of the frictional force is proportional to the mean value of —Svf \ at low impact 
parameters, the distribution of Svf x 's is shifted toward zero and the frictional force is reduced. 




Fig. A7. — Dynamical friction reduction factor as a function of angle between the normal to the binary's orbital plane and the relative 
velocity vector of the field stars, for V l~Vun = 1- Different curves correspond to different values of the maximum impact parameter; proceeding 
upward, p m ax/a = (0.2, 0.5, 1, 2, 3, 5, 20). Points are averages computed from the numerical integrations; curves are spline fits. 



24 



Brownian Motion 




Merritt 



25 




Fig. A9. — Dimensionless ratio D3 (V, p m ax) (eq. 
to the initial velocity, normalized to the value for 



point-mass scatterer 
lowest (highest) V produces the highest (lowest) D3 at small p m ax- 



72t ) describing changes in the mean square field-star velocity in a direction perpendicular 
Different curves correspond to V/Vti n 



0.1,0.2,0.5,1,2,3; the 
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r max 



Fig. A10. — Increase in the diffusion coefficient (Ai> 2 ) for a massive binary moving at low velocity, V-^~ °7> compared to its value for a 
point mass. Curves are labelled by the value of Vti n /uf. Dashed line is the approximation of equation (J74[) . 



Merritt 




